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ABSTRACT 

Aims. We show that observation of the time-dependent effect of microlensing of relativistically broadened emission lines 
(such as e.g. the Fe Ka line in X-rays) in strongly lensed quasars could provide data on celestial mechanics of circular 
orbits in the direct vicinity of the horizon of supermassive black holes. 

Methods. This information can be extracted from the observation of evolution of red / blue edge of the magnified line 
just before and just after the period of crossing of the innermost stable circular orbit by the microlensing caustic. The 
functional form of this evolution is insensitive to numerous astrophysical parameters of the accreting black hole and of 
the microlensing caustics network system (as opposed to the evolution the full line spectrum). 

Results. Measurement of the temporal evolution of the red / blue edge could provide a precision measurement of the 
radial dependence of the gravitational redshift and of velocity of the circular orbits, down to the innermost stable circular 
orbit. These measurements could be used to discriminate between the General Relativity and alternative models of the 
relativistic gravity in which the dynamics of photons and massive bodies orbiting the gravitating centre is different 
from that of the geodesics in the Schwarzschild or Kerr space-times. 

Key words. Gravitational lensing: strong; Gravitational lensing: micro; quasars: supermassive black holes; Black hole 
physics 


1. Introduction 

General Relativity (GR) is by far the most popular theory 
of the relativistic grav ity, tested w ith high precision via a 
set of “classical” tests ( [Will |2Q14 ) in its weak field regime, 
when the gravitational potential is 
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where M and R are the mass and size of the gravitating 
system (we use the system of units in which the speed of 
light c = 1). However, none of the classical tests of GR pro¬ 
vides a verification of the theory in the strong-field regime 
U- 1. 

At the same time, several observational phenomena oc¬ 
cur in this regime. The U ^ 1 regime of gravity is infiuenc- 
ing the physical processes in the vicinity of neutron stars 
and stellar mass black holes in pulsars, in the X-ray binaries 
and in the cores of massive stars exploding as supernovae. 
Electromagnetic emission from supermassive black holes in 
the centre of the Milky Way and other galaxies is also af¬ 
fected by the strong gravity effects. Einally, the dynamics of 
expansion of the Universe within the Hubble distance Rh 
occurs in the strong field regime with U ^ GnPoR^h ^ ^ 
(here po — 10“^^ g/cm^ is the critical density of the 
Universe). 

Although the relativistic gravity effects certainly influ¬ 
ence the observational proper ties of the sources powered by 


ence the observational proper ties oi the sources p < 
neutron stars and black holes (|Eabian et al. 2QQQ| Remillard 
fc McClintock | [2QQ^ Psaltis ' 2QQ8p , large uncertainty of 


the physical conditions and parameters describing physical 
phenomena in these systems does not allow to test the GR. 
The unique nature of the expanding Universe system also 
precludes the possibility to test the theory. The dynamical 
models of the Universe do not agree with observations in a 
straightforward way, as it is clear from the necessity of in¬ 
troduction of the Dark Matter and Dark E nergy in certain 
fine-tu ned proportions 0.3, ^ 0.7 (Planck Collab. 

[2015 ) into the description of the dynamics of the Universe. 


GR is not the uniqu e viable theore t ical model of rela¬ 
ti vistic gravity (see e.g. |Clifton et al. (|2012P; |Berti et al. 


(2015) for reviews). Some of the alternative theories are 


motivated by the possibility to solve the problem of the 
Dark Matter and Dark Energy. Others are introduced to 
achieve compatibility of the relativistic gravity with quan¬ 
tum mechanics. All the extensions and modifications of GR 
have the same weak field limit as the GR and as such are 
compatible with the classical tests of GR. At the same time, 
predictions of the extensions / modifications of GR in the 
strong field regime could strongly deviate from the predic¬ 
tions of the GR. In particular, the dynamics of the Universe 
and the properties of the black holes in different alternative 
models of relativistic gravity could be significantly differ¬ 
ent. 

Perhaps, the most straightforward possibility to test dif¬ 
ferent properties of black holes in the alternative models of 
relativistic gravity is via a measurement of the laws of ce¬ 
lestial mechanics in the strong gravitational field. After all, 
it is the study of the celestial mechanics of planets in the 
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Fig. 1. Map of the red/blue shift of an emission line from 
the equatorial disk around black hole with a = 0.9 inclined 
at 80° with respect to the line of sight. Green thin contours 
show the radial coordinate of the point of emission from the 
black hole in the Boyer -Lindquist coordinates. Contours 
start from r = 2 with the increment of Ar = 2 up to r = 
20. Cyan thick contours show the redshift factor, starting 
from 0.3 with increments of 0.1. The trajectory of caustics 
inclined with respect to the line of nodes of the disk is 
shown by the red line. 


Solar system that led to the formulation and verification of 
the non-relativistic theory of gravity in the past. 

Several possibilities exist for such measurements using 
the observations of the supermassive black hole in the cen¬ 
tre of the Milky Way. Observations in the millimetre wave¬ 
length range have reached the angular resolution of about 
10 /ias necessary to resolve the event horizon of the su¬ 
permassive blac k hole in the Galactic Centr e and several 
nearby galaxies ( Doeleman et al. [2008 


2012). Precision of 


the astrometric measurements with the infra-red interfer¬ 
ometric observations will also reach a compara ble angular 
scale with the GRAVITY instrument at VLT dCliessen et 

HTIPMol ). 


The angular size scales of more distant supermassive 
black holes in the galaxies beyond the Milky Way are typ¬ 
ically too small to be directly resolvable. For these distant 
black holes a potential possibility for the study of strong 
gravity effects is via the measurement of the Doppler and 
gravitational redshift of atomic lines emitted from the in - 
nermost portion of the accretion disks (Fabian et al. 120^ 


However, if such measurements are not complemented by 
the information on spatial location of the portions of the 
disk with certain red / blue shifts, no real study of “celestial 
mechanics” of the orbits around the black hole is possible. 

In what follows we propose a method for the spatial lo¬ 
calization of the red / blue shifted parts of the accretion 
disks. The method is based on the observation of the effect 


ally lensed quasars ( 

iDai et al. 2010 |Ghen et al. 2012 

Ghartas et al. ||2012 

|MacLeod et al.||2015|). The microlens- 

ing caustics, passing in front of the accretion disk, work as 


a “scanning magnifying glass”, helping to resolve the region 
of emission of the relativistically broadened lines. We argue 
that measurement of the time evolution of the red / blue 


edges of the line, magnified by the caustics, provides the in¬ 
formation on the radial dependence of the rotation velocity 
and gravitational redshift down to the location of the in¬ 
nermost stable circular orbit (ISCO) around the black hole. 
This opens a possibility to test of the GR in the strong field 
regime. 


2. Relativistically broadened lines from the black 
hole accretion disks 

Radiatively efficient accretion onto supermassive black 
holes in quasars and Seyfert galaxies is conventionally as¬ 
sumed to proceed through an accretion disk which rotates 
at nearly Ke plerian velocities in the equatorial plane around 


a black hole ( Shaknra fc Sunyaev |1973" Novikov & Thorne 


1974). Photons emitted by the matter moving along the 
circular orbits are gravitationally redshifted and also red 
/ blue shifted by the Doppler effect. The red and blue 
shifts depend on the position of the emission point and 
on the direction of photon propagation. Superposition of 
the red/blue shifted atomic line emission from different po¬ 
sitions in the accretion disk leads to t he occurrence of rel¬ 
ativistically broade ned emission lines ( Tanaka et al. 

Fabian et al. [2000 ). 


Galculation of the spectral shape of the broad line can 
be done following the standard procedure o f “ray shooting” 
from the observer towards the black hole (Bromley et al. 


1. |1997D . Fig. shows an example of such 
particular case of the black hole with the 


1996; Fanton et al. 
calculation for a pa'rticui 
specific angular momentum a = ^.^GnM (M is the black 
hole mass). The image is calculated by tracing the trajecto¬ 
ries of photons emitted along the “line of sight”, inclined at 
i = 80° with respect to the black hole rotation axis. Photon 
trajectories start at a remote “screen”, placed perpendicu¬ 
larly to the line of sight. Photons are back-traced from the 
“screen” till the moment when their trajectories intersect 
the equatorial plane, which contains the accretion disk, or 
when they hit the black hole horizon. Photons crossing the 
equatorial plane M times form the Afth. order image (only 
the first order image is shown in Fig. f^. Photon trajecto¬ 
ries are calculated using the Runge-KS;ta method for the 
solution of the first-order differential equations of the null 


geodesics (Bardeen et al. 1972). Only photons which hit 


the equatorial plane at the distances larger than the dis¬ 


tance of the ISGO ( Bardeen et al. |[l972 ) are contributing 
to the image. 

The rotation velocity of the prograde circular orbit at 
the distance r in the equatorial plane around the black hole 
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in the Locally Non-Rotating Reference Frame (LNRF), de¬ 
fined by the tetrad 
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A sin^ 0, = r‘^ + cos^ 0 and 


(Bardeen et al. Il972| ). 


where A 
A = - 2R 

The Doppler factor of the photon emitted at an angle 
aem with respect to the direction of the orbital motion in 
the LNRF is 


l-VWcOs{aem) 

Energy of the photons measured by an observer at infinity, 
Eobs^ is related to the reference line energy Eem through 
the conserved zero component of the photon momentum 
p^. It is expressed in the Boyer-Lindquist coordinates as 


Eobs — Pt — ^P{a) 


( 5 ) 
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where is the photon four-momentum in the LNRF. 

In the image of the strongly inclined accretion disk in 
the left panel of Fig. one can clearly see the Doppler effect 
of the disk rotation, which leads to the red shift of the right 
part of the disk image, where aem ~ tt, and to the blue 
shift of the left part of the image, where aem ~ 0- From 
Eq. <§ one can see that the red / blue shift dependence on 
the angle of emission is through the conventional Doppler 
factor in the LNRF and additionally through the effect of 
the LNRF dragging by the black hole rotation. 

The energy profile of the relativistically broadened line 
can be obtained via a convolution of the maps of 5obs = 
(Eobs/Eem) (see Fig. [B with the radial intensity profile of 
the line emission /em(n- 


E{Eobs) = / / < 

J Ja,^ 




obs 
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where a^ (3 are the coordinates on the observer’s “screen”. 
Examples of such convolution for different radial pro¬ 


files le 


and different outer radii of the line emis¬ 


sion region are shown in Fig.|^ One can see that the func¬ 
tional shape of the broad line spectrum depends on the 
parameters like the radial brightness profile and the size of 
the line emission region. These parameters are not directly 
measured. Instead, they have to be derived from the obser¬ 
vational data. It is clear that model fits to a rather wide 
range of spectral shapes could be found by the appropriate 
adjustments of the functional parameter 

Thus, even though the spectral shape of the relativisti¬ 
cally broadened line is sensitive to the details of the rela¬ 
tivistic gravity theory in the strong field regime, measure¬ 
ments of the spectral shape alone can not be used to test 
the relativistic gravity theory. One of the discussed possi¬ 
bilities for the tests of the strong gravity using the data 
on the relativistically broadened lines is to complement the 
line spectrum measurements with the additional timing in¬ 
formation, available e.g. in th e out-of-steady-state config¬ 


uration of the accretion disk (Bambi 2013 Johannsen & 


Psaltis 2013). In the following section we propose an al 


ternative type of timing analysis, applicable in the case of 
the broad line emitting system which is a part of the strong 
gravitational lens. 



Fig. 2. Spectra of the broad emission line from a region 
of the size Rdisk with radial brightness profile r~^ in the 
disk around a black hole with specific angular momentum 
a = 0.9, inclined at i = 45° with respect to the line of sight. 


3. Microlensing of the relativistically broadened 
lines in strongly lensed quasars 

3.1. Qualitative picture 


Quasars emitting relativistically broadened emission lines 
can occa sionally be parts of the strong gravitational lens 
systems (Schneider et al. |1992 Blandford & Narayan 1992| 


Claeskens fc Snrtej |2002p ,~ in which case multiple images of 

the quasar may appear. Typical deflection of light rays by 
the leasing galaxy of the mass Ml is by an angle 
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In addition to the strong macro-lensing effect, the source 
images are affected by the effect of mic rolensing by mas- 


sive sub-structures 

in the lensing galaxy (Chang & Refsdal 

19791 IKayser et al. ||19861 |Kochanek |2004| Schmidt & 

Wambsganss ||2010 

). The most commonly considered type 


of microlensing is by stars, i.e. point-like objects with the 
masses m ^ Mq. An individual star significantly magnifies 
the flux from the background sources situated within its 
Einstein ring 
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Here dL^ds^dLs are the angular diameter distances to the 
lens, to the source and between the lens and the source. 
d ^ dL^ds^dLS is the distance scale. The projected linear 
size of the Einstein ring on the source plane is 
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Fig. 3. Example of the network of the microlensing caustics 
for the microlensing optical depth r 0.7 for the region 
approximately 2.5x 10^^ cm by 1.5x10^^ cm in size. Colour 
scale represents the magnification. White box marks the 
region of caustics which we used for the calculations of the 
effect of caustic crossing on the line emission. 


The microlensing selectively magnifies the flux from sub¬ 
structures with the linear size Rs smaller than r^;. In this 
case the magnification factor can be estimated as 


/^mZ 


Rs 


1/2 


( 10 ) 


If the lens galaxy mass in its innermost part within 

Re = dL^E ^ 6 [Ml/IO^^Mq][ d^/l Gpc]^^^ kpc is 
dominated by stars, the total number of the micro lenses 
inside the macro-lensing Einstein ring of the radius &e is 
Nmi ^ Ml/tti. The Einstein circles of individual stars over¬ 
lap in projection on the sky, because the total area covered 
by the Einstein rings of the individual stars is comparable 
to the area of the macro-lensing Einstein ring: 


'^ml — 


Nmiirdl 
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7r0| 


~ 1 


( 11 ) 


The large, Tmi ^ 1, “microlensing depth” leads to a com¬ 
plicated pattern of (de) magnification and to formation of 
a network of caustics. An example of such a network is 
sh own in Eig. calcu l ated using the ray-t racing method 


of Wambsganss (1990); Bate et al. ( 2Q1Q[) for Tm i =0.7 


and the initial stellar niass function oT Chabrier | ( 2003[ ). 
The magnification factor of a point source situated exactly 
behind the caustic curve is infinite. The maximal magnifi¬ 
cation factor for a source of finite size Rs can be estimated 
from Eq. 

The microlensing magnification factor is variable in 
time, because of the relative motion of the source, of 
the lens, of stars in the lens galaxy and of the observer. 
Typical estimates of the relative projected velocity are 
V ^ 10^ km/s. A source of the size Rs passing behind the 
caustics stays magnified during a period of time 



\ 1 
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3 X 10^"^ cm 

10^ km/s 


n -1 


month (12) 


The time variability effect of the microlensing is now 
observed in a range of the strongly lensed quasars, via 


the measurement of the difference in th e lightcurves of the 
multiple images of the lensed sources ( [Chang fc RefsdS 
||1979[ [Schmidt fc: Wambsganss ||2Q1Q[ [Dai et al. ||2Q1Q| 

^ "2Q12| [Chartas et al. [[2Q12| [MacLeod et al!] 


Chen et al. 


2015). Measurements of the strength of the microlensing 


visible ([Kochanek [[2004 Eigenbrod et al. 2008), and X- 

ray (Dai et al. 120101 

Chen et al. |2012| Chartas et al. |2012 

MacLeod et al.[[2015 

) emission in the lensed quasars. 


3.2. “Microlensing portrait” of the relativistically broadened 
line 


If accretion disk is a source of the relativistically broad¬ 
ened line emission, the line flux is also subject to the mi¬ 
crolensing effect. Strong magnification of the line flux is ex¬ 
pected during the periods of caustics crossings of the disk 
images ( [Popovic et ah, [[2006 Chartas et al. 2012). 

In general, the caustics crossing does not affect equally 
the flux of the “red” and “blue” wing of the line. Erom 
Fig.g one could guess that the caustic, passing in front of 
the disk image, say from right to left would first magnify 
the flux of the red wing of the line and then the flux of the 
blue wing of the line (the only exception is when the caustic 
is parallel to the line of nodes of the disk). This expectation 
is confirmed with the numerical simulation shown in Eig.|^ 
The left panel of the figure shows an accretion disk inclined 
at i = 45° and truncated at the radius r = 20GArM. The 
disk surface brightness profile is assumed to be /(r) (x 
The right panel of the figure shows the trajectory of the 
disk passing behind a caustic taken from the map of Eig.[^ 
The motion of the disk is along the dashed line, so that the 
caustic passes perpendicularly to the line of nodes of the 
disk. 

The middle panel of Eig. [^ shows the evolution of the 
spectrum of the relativistically broadened line during the 
period of caustic crossing. The y axis shows the energy in 
the units of the reference line energy. The x axis shows the 
position of the caustic centre with respect to the black hole 
and effectively represents the time. The colour scale shows 
the intensity of emission at different energies and different 
moments of time. One can see that the amplification of 
the flux from the part of the disk which is exactly behind 
the caustic leads to the appearance of an increased flux 
between some minimal and maximal energy. The increased 
flux component has pronounced red and blue “edges”. The 
energy of the “edge” is determined by the energy of photons 
^obs (Efl- E’ calculated for the photons emitted from the 


— point of intersection of the caustic with the line of nodes 
of the disk, if the caustic does not pass through the 
ISCO or 

— points of intersection of the caustic with the ISCO, if 
the caustic is in front of the ISCO. 


The energy of the edge evolves with time because the points 
of intersection of the caustic with the line of nodes of the 
disk and with the ISCO move as the caustic advances across 
the disk image. 

In principle all the Afth. order lensed images of the disk 
contribute to the line flux. However, the flux from higher 
order images is typically weaker than that of the first order 
image. The bottom panel of Eig. [^ shows the evolution of 
the line spectrum produced by the microlensing of the sec¬ 
ond order image. One could see that the evolution of the 
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Fig. 4. Top left: image of the red / blue shifted emission from an accretion disk with external radius 20Rg around a 
black hole with specific rotation moment a = 0.9Rg inclined at i = 45° with respect to the line of sight. Color scale 
gives the red / blue shift Eobs/Eem- Top right: a microlensing caustic and its trajectory (dashed line) in front of the 
accretion disk image (shown by a solid contour). Colour scale represents the log^g of the microlensing magnification. Top 
middle: evolution of the spectrum of relativistically broadened line during the period of caustic crossing in front of the 
disk image, x axis is the distance of the point of intersection of caustic with the line of nodes of the disk, measured in 
the units of Rg. Only first order image is taken into account. Color scale gives the relative intensity in arbitrary units. 
Bottom: the same as the top images, but for the second order image of the disk. 


second order spectrum is significantly different from that of 
the first order image. If the sensitivity of telescope is suf¬ 
ficient, the spectral evolution should reveal multiple edges 
corresponding to the images of different order. 

The shape of the curve which the red and blue edge 
draw in the middle panel of Fig. depends on two pa¬ 
rameters: the rotation moment of the black hole a and the 
inclination angle of the disk i. The shape does not depend 
on numerous additional parameters of the physical system 
harbouring the black hole, the accretion disk and the X-ray 
source illuminating the accretion disk. It also does not de¬ 
pend on the parameters of the microlensing caustics, such 
as the projected velocity of the caustic, the angle of incli¬ 
nation of the caustics with respect to the line of nodes of 
the disk, the typical mass of the stars in the lensing galaxy 
and their distribution in masses etc. It is rather the time 
dependence of the full spectrum of the line which is sensi¬ 
tive to the numerous astrophysical parameters. A range of 
the edge time evolution curves for different values of a and 
i parameters is shown in Fig.l^ In a sense, these red / blue 
edge time evolution curves could be called the “microlens¬ 
ing portraits” of the broad emission lines. 

The independence of the shape of the microlensing por¬ 
trait of the line on the system parameters could be un¬ 
derstood in simple terms if one uses a toy model of the 


accretion disk based on the Newtonian gravity instead of 
the GR (or an alternative theory of gravity). This toy 
model neglects the gravitational redshift of photons, taking 
into account only the Doppler effect. The circular orbits in 
Newtonian gravity have velocities 


VK = ^ (13) 

The maximal Doppler red / blue shift of line emission pro¬ 
duced at the distance r occurs for the emission produced 
at the line of nodes of the disk. It is determined by the 
inclination angle i of the disk 


/ N /-N I GnM cos‘^ i . 

g{r) = vk cos(2) = Y- - - (14) 

The microlensing caustic is a line inclined at an angle ac 
with respect to the line of nodes of the disk. If the caustic 
moves with the projected velocity Uc, the motion of the 
intersection of the caustic with the line of nodes is described 
by the relation 

r(t) = Uct cos Q^c (15) 
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Fig. 5. Evolution of the energy of the red / blue edges of 
the line throughout the caustic crossing period for different 
rotation momenta of the black hole and different viewing 
angles. 


so that the time evolution of the maximal / minimal blue 
/ red shift is 


9{t) = ± 


GatMcos^ i 

Vc COS a 


(16) 


The functional shape of g{t) is 1/Vi and it does not depend 
on the parameters of the system. In fact, in non-relativistic 
theory the functional shape of the curve is determined by 
the gravity theory itself, rather than by the parameters of 
the system. If the gravity force law would deviate from the 
1/r^, the functional shape of the g{t) curve would deviate 
from l/y/t. 

Variations of parameters of the system do not change 
the shape of the curve, but instead they influence the scal¬ 
ing of the X and y axis in the “microlensing portrait” of the 
curve. For example, changing the velocity of motion of the 
caustic would lead to stretching or squeezing of the time 
scale of the caustic crossing, i.e. the scale of the x axis. 
A similar effect occurs if the angle ac changes. Only the 
projected velocity of the caustic, Vc cos a, matters. Increase 
or decrease of the mass of the black hole also changes the 
scaling of the x axis, because the natural scale of variations 
of the red / blue shift is given by the gravitational radius 
Rg = G ]\f M. 

Description of the “microlensing portrait” of the line 
within GR introduces only technical complications, but 
does not change this qualitative picture. The red/blue shift 
at the point of intersection of caustic with the line of 
nodes or with the ISCO can be calculated from Eq. (§, 

9 — ^obs / ^em 1 - 

The parameters r, aem in @ , which are functions of 
time, have to be calculated numerically via back-tracing of 
photons from infinity to their emission points in the disk. 


The asymptotic form of these functions for the times when 
the intersection of the caustic with the line of nodes of the 
disk is far away from the black hole is given by aem — 7r/2—i 
and by Eq. (15). 

The redshift g has two contributions. In addition 
to the Doppler effect given by the Doppler factor 5emi 
the line emission is gravitationally redshifted. The func¬ 
tional shapes of both the gravitational redshift and of the 
Keplerian velocity depend on the black hole rotation mo¬ 
mentum a. 

One more important difference between non-relativistic 
and relativistic gravity model is that in non-relativistic 
gravity the function g(t) is not bounded from above and 
from below. This is because the circular orbits can be arbi¬ 
trarily close to the centre of gravity. To the contrary, in GR 
the stable circular orbits exist only down to the distance of 
the ISGO. Finite gravitational and Doppler red/blue shifts 
of photons emitted from the ISGO explain the “plateau” 
parts of the “microlensing portraits” of the line shown in 
Fig.|5] 


3.3. Celestial mechanics of circular orbits around black hole 
from the microlensing portraits of emission lines 

Fig. 1^ shows that the edge evolution curves (the microlens¬ 
ing portraits) form a family of curves parametrised with the 
black hole rotation moment a and the disk inclination angle 
i. The scaling of the x and y axes depend on two additional 
parameters: the angle between the caustic and the line of 
nodes of the disk a and the velocity u, which provides the 
relation x = vt between the time t and the x coordinate in 
Fig-i Overall, the full curves depend on four parameters. 
These four-parameter curves could be fitted to the obser¬ 
vational data, presented in the same form as the simulated 
data shown in the middle panel of the Fig. Note that in 
spite of the fact that the overall evolution of the form of the 
spectrum of the broad line throughout the caustic crossing 
period is quite complicated, the evolution of the edge could 
be determined rather reliably and independently from the 
details of the complex evolution of the spectrum. Taking 
into account sharpness of the edge, one could use the ded¬ 
icated imaging analysis methods (e.g. wavelet analysis) to 
single out the edge evolution on top of the overall spectral 
evolution. 

If the quality of the data is high enough, fitting of the 
model curves to the data provides a test of the relativistic 
gravity theory (if the number of data points in the edge 
evolution curve is larger than four). If the sensitivity of 
the telescope is high enough, even more constraints on the 
gravity theory could be derived from the observation of the 
edge evolution curves of the higher order images (see the 
bottom panel of Fig. [^, because the shape of these curves 
does not depend on any additional parameters. 

Depending on the orientation of the accretion disk with 
respect to the line of sight, the edge evolution curves pro¬ 
vide different types of information on the details of the 
structure of gravitational field and of the motion of mas¬ 
sive test bodies and photons in this gravitational field. 

If the disk is observed almost face-on (the upper panel 
of Fig. [^, the edge evolution curve has two red branches 
{EobslEem < I) and no blue branch. This is obviously due 
to the effect of the gravitational redshift. The asymptotic 
behaviour of the two red branches at large x = vt is almost 
identical. Taking the moment of time with respect to which 
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the two edges are symmetric, one could define the moment 
when the caustics crosses the centre of the disk image. If the 
time is counted with this moment taken as a zero reference 
time, the shape of the g{t) curve provides a measurement 
of the gravitational redshift as a function of the distance 
from the black hole. 

The “plateau” of the edge evolution curve, visible in the 
upper panel of Fig.[^ corresponds to the time interval when 
the caustic passes the ISCO. The depth of the plateau is 
determined by the gravitational redshift of photons emitted 
almost perpendicularly to the disk from the ISCO. This 
redshift depends on the black hole rotation moment a. The 
width of the plateau is determined by the size of the ISCO 
which is also a function of a. 

If the disk is significantly inclined, the Doppler effect 
“tilts” the microlensing portrait, lifting one side of the edge 
curve up and pushing the other side down energy-wise. 
Obviously, the measurement of the distance (time) depen¬ 
dent tilt of the curve provides a measurement of velocities 
of the circular orbits as a function of the distance. 

If the black hole rotation moment a is not zero, the incli¬ 
nation of the disk not only “tilts” the plateau correspond¬ 
ing to the period of passage of the ISCO, it also deforms its 
shape. This is because of the fact that even in the absence 
of the Doppler effect due to the disk rotation, there is an 
angular dependence of the gravitational redshift in Eq. (§ 
caused by the rotational drag by the black hole. Photons 
reaching the observer from different points of the ISCO 
are emitted in different directions spanning a cone with the 
axis given by the black hole rotation axis and with an open¬ 
ing angle comparable to the inclination angle of the disk. 
Evolution of the energy of the red / blue edge of the line 
during the period of caustic crossing of the ISCO provides 
a “scan” of the red / blue shift of the photon energy as a 
function of direction along the cone. 

If the quality of the spectral measurements is high 
enough, the angle between the caustics and the line of nodes 
of the disk could be measured from an effect of “splitting” of 
the edge during the period of caustic passage of the ISCO. 
The origin of such splitting is clear from Eig. The red 
/ blue shifts of photons emitted from the points of inter¬ 
section of the caustic with the ISCO are different. If the 
caustic is perpendicular to the line of nodes, this difference 
is small and the effect of the splitting of the edge is difficult 
to observe. However, if the caustics is strongly inclined, the 
effect of the splitting by the caustic is more pronounced 
and could be used for the measurement of the angle a. 

4. Discussion 

The results of the previous section show that measurements 
of the evolution of the red / blue edge of line emission from 
the accretion disk during the periods of crossing of the disk 
image by the microlensing caustic provide a unique possibil¬ 
ity to make a range of measurements, which could be used 
to characterise the structure of the strong gravitational field 
in the direct vicinity of the black hole horizon. 

The functional shape of the microlensing portraits of the 
lines, shown in Eig. is largely independent from the pa¬ 
rameters of the accretion disk and of the caustic passing in 
front of the disk image. This near-parameter-independence 
indicates that the measurements of the “microlensing por¬ 
traits” of the lines could potentially be used to test the 
GR in the strong field regime near the black hole horizon. 


Measurements of the “microlensing portraits” would, in a 
sense, be equivalent to the measurements leading to the 
Kepler laws, on the basis of which the celestial mechanics 
of the Solar system was established and the non-relativistic 
gravity law was derived in the past. 

The most straightforward observational example of the 
relativistically broaden ed line is the Ee Ka line at the ref- 

Broad Ee Ka 


2000 ). 


erence energy 6.4 keV (Eabian et al. 
line is observed in a large number of Seyfert galaxies and 
quasars, starting from the initial A SCA discovery of the line 
in MCG 15-6-30 galaxy spectrum ( Tanaka et al. | 1995). 

A number of strongly lensed quasars, such as e.g. 
RX J113 1-1231, possess the broad Ee Ka line in thei r 


spectra (Ghartas et al. 2012 Reis et al. 2014). 


Moreover, Ghartas et al. ( 2012| ) find that the broad Ee Ka 
line emission of RX J1131-1231 is influenced by the mi¬ 
crolensing. They interpret the shifts in the energy of the 
peak of the line as the magnification due to the microlens¬ 
ing effect. Unfortunately, th e time coverage of the observa¬ 
tional campaign reported in [Ghartas et al. (2012) was not 
sufficient for the measurement of the evolution of the line 
(or edge) energy with time. An estimate of the black hole 
mass in RX J1131-1231 is M - lO^M©, so that the caustic 
crossing takes approximately one month (see Eq. 12). This 
time was also the typical time interval between the sub¬ 
sequent observatio ns in the observationa l campaign of the 
source reported by Ghartas et al. (2012). Only two obser¬ 
vations are available for the entire caustic crossing period. 

In spite of the fact that the observations of RX J1131- 
1231 do not allow the measurement of the “microlensing 
portrait” of the Ee Ka line, they are useful for the demon¬ 
stration of the feasibility of such measurement. The shift 
of the line energy due to the microlensing was observed in 
the weakest component “D” of the source image. Brighter 
images are also known to be affected by the microlensing 
variability (|Dai et al. 2010). Monitoring of the source in 
X-rays should sooner or later catch the caustics crossing 
periods in the brighter components of the source image. 

The caustics crossings of the image components of the 
RX J1131-1231 and of other strongly lensed quasars should 
be frequent enough (an estimate of “s everal” per decad e 
could be derived for RX J1131-1231 from Dai et al. (2010)). 
An efficient observational approach for the measurement of 
the “microlensing portrait” of the Ee Ka line should in¬ 
clude also an efficient “forecasting” of the caustics crossing 
periods, followed by the intensive observational campaign 
with the dense time coverage in X-rays, upon the positive 
outcome of the forecast of the caustic crossing period. 

The forecasting of caustics crossings could be done using 
the source monitoring campaigns in X-rays or in the visible 
band. The overall size of the X-ray emitting part of the 
disk is just about 10 gravitational radii. It is typically much 
smaller than the Einstein radius of the microlensing centres. 
This results in strong magnification of the X-ray ffux of 
individual images of the source, by up to one-two orders of 
magnitude during the periods of the caustics crossings. 

Such strong increase of the flux of one of the source 
images could be detected via a monitoring campaign of a 
strongly lensed source with an X-ray telescope with even 
a moderate angular resolution, which is not able to resolve 
the multiple images of the source (like e.g. the Swift/XRT 
instrument). Erom the pattern of the magnification maps in 
Eig. one can see that the caustic crossing moment should 
be preceded by a period of the steep rise of the magni- 
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fication factor in about 50% of caustics crossing episodes 
due to the presence of the “plateau” regions in the caus¬ 
tic shapes. Identification of such steep flux rises in caustic 
crossing precursors with a regular monitoring campaign of 
a lensed source (year-to-multi-year long) should pose no 
problem. 

The size of the accretion disk in the visible band is sig¬ 
nificantly larger tha n the gravitational radius (by one-two 
orders of magnitude, [Dai et al. |2QlQ[[Eigenbrod et al. |2QQ8| 


MacLeod et al. 2015). The overall changes of magnitudes 
of the source iniages due to the microlensing are moderate, 
because of the large source size. Still, the caustic crossing 
periods should be identifiable as the periods of the rapidly 
growing amplitude of the difference of the flux of one of the 
source images. The best possibility for the caustics cross¬ 
ing detection is, perhaps, a combination of the data of the 
visible and X-ray band monitoring. Optical monitoring of 
the source could serve for the confirmation of the approach 
of the caustic crossing moment, identified based on the de¬ 
tection of the “X-ray precursor”. 

As soon as the caustic crossing precursor is detected, an 
intense observational campaign with more powerful X-ray 
telescopes which are able to resolve the individual source 
image components could be initiated, with the main goal to 
obtain a high signal-to-noise ratio “microlensing portrait” 
of the relativistically broadened line, like the one shown in 
the middle panel of Fig. 

The quality of the test of GR, which could be obtained 
via the measurement of the microlensing portrait of the 
Fe K(a, line depends on the precision of the measurement 
of the edge evolution curve. This precision depends on the 
number of the observations, carried out throughout the 
caustic crossing period, and on the precision of the measure¬ 
ment of the edge position in each observation. Minimum five 
observations per caustic crossing are required to fully de¬ 
termine the parameters of the curve. Each additional data 
point is potentially useful for the test of the theory. 

First constraints on the relativistic gravity theory could 
be obtained using the proposed method already with the 
current generation X-ray telescopes - Swift/XRT, XMM- 
Newton, Chandra and Suzaku. Higher precision of the test 
can be obtained with the high spectral resoluti on measure- 
ment, e.g. the X-ray spectrometers on Astro-H dMitsuda et 
id~||2Qi4l ) and/or Athena ( |den Herder et alT ]|20Wf^ 
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